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1. Introduction 

A further dimension to high time resolution astrophysics (HTRA), largely unexplored so far 
[38, 31], is rapid changes in the optical polarization of astronomical sources. However, with new 
[26, 30, 8] and future instrumentation, this new region of parameter space may provide vital clues 
to aid us in understanding the physics of some variable objects. Here, we focus on variability of 
polarization in low-mass X-ray binaries (LMXBs), and show that rapid variability is present and 
can be detected in the brightest sources with current instrumentation on moderate size telescopes. 
A global picture of polarization and magnetic properties of the accretion flows in these objects 
will come essentially via studies of populations of LMXBs using future facilities that will provide 
improvements in both sensitivity and temporal domain. 

Optical and infrared light from LMXBs can appear polarized either due to intrinsic polarization 
of the emitting photons, or by the (Thompson or possibly Rayleigh) scattering or absorption of 
unpolarized photons. Polarization of optical light from LMXBs due to absorption by interstellar 
dust is inferred by a characteristic [37], unchanging, wavelength-dependent polarization level [39]. 
Scattering of intrinsically unpolarized thermal emission by electrons [11] within the system has 
been detected as optical polarization that changes as a function of the orbital phase of the binary in 
two black hole (BH) LMXBs, A0620-00 and GRO J1655-40 [12, 18]. The reported strength of the 
polarization from scattering is typically a few per cent. All other reported measurements of optical 
polarization in LMXBs have an interstellar dust origin as far as we are aware, except possibly for 
the BH XTE Jl 1 18+480 [39] and the neutron star (NS) LMXB Aql X-l [7]. 

Very recently, evidence for polarization intrinsic to the emitting photons of LMXBs has been 
revealed from infrared observations [40, 32]. Shahbaz et al. (2008) [40] performed infrared spec- 
tropolarimetry of three LMXBs and found two of them (the NSs Sco X-l and Cyg X-2) to be in- 
trinsically polarized, with an increasing linear polarization (LP) at lower frequencies (up to ~ 2% 
in Sco X-l and ~ 10% in Cyg X-2). This dependence on wavelength cannot be explained by 
scattering or interstellar dust absorption. 

We performed near-infrared (NIR) imaging polarimetry of six LMXBs in 2005-6 using the 
Infrared Spectrometer And Array Camera (ISAAC) on the Very Large Telescope (VLT) (European 
Southern Observatory [ESO] programme 275.D-5062) and the UKIRT 1-5 micron Imager Spec- 
trometer (UIST) + Infra- Red Polarimetry facility (IRPOL2) on United Kingdom Infrared Telescope 
(UKIRT) (making use of the UKIRT Service Observing Programme; UKIRTSERV). The observa- 
tions and results are summarized in Table 1; this work was published in Russell & Fender (2008) 
[32]. Most sources were in quiescence, and no polarization was detected. LP was detected (at the 
> 3a level) in two active sources, the NS Sco X-l and the BH GRO J1655-40 (Fig. 1). The LP 
is strongest in the reddest NIR filters, and is variable for Sco X-l. For this latter source, data were 
acquired on four dates and the LP was seen to vary both between dates and within the same night 
(the time resolution was 2.5 minutes). 
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Table 1: Details of the 2005-6 NIR polarimetric observations of LMXBs taken with the VLT + ISAAC 
(GRO Jl 655-40) and the UKIRT + UIST (the other five), and main results, lcr errors are given for LP, or 
3<7 upper limits. 



Source 


Epoch (days) 


Source activity 


Filters 


LP (filter) 


GRO J0422+32 


2006 Aug (1) 


Quiescence 


J,K 


< 12% (J) 


4U 0614+09 


2006 Oct (1) 


Persistent 


K 


< 16% (K) 


XTEJ1 118+480 


2006 Feb (2) 


Quiescence 


J,K 


< 6% (J), < 7% (K) 


Sco X-l 


2006 Feb-Mar (4) 


Persistent 


J,H,K 


0.29-0.65 (J), 0.07-0.57 










(H), 0.13-0.91 (K) 


GRO J1655-40 


2005 Oct (2) 


Hard state 


J,H,K S 


4.7±1.5 (H), 5.8±1.9 (K s ) 


Aql X-l 


2006 Aug-Oct (3) 


Quiescence 


J,H,K 


< 5% (J), < 17% (H), 










< 12% (K) 




Figure 1: NIR intrinsic LP detected in Sco X-l and GRO J1655-40, as a function of frequency. For 
Sco X-l we include the optical data (filled circles) and the fit to these data for an interstellar origin to the 
polarization (solid curve) from [39]. The crosses on the left with vertical bars are our KHJ-band results 
(2(7 detections of LP); a clear variable LP component above that expected from interstellar dust is seen 
at the lowest frequencies. For GRO J1655-40, VRI-band data are from [18] and a K-band upper limit in 
quiescence is from [13]. The lines in the lower right panel show the fit to the outburst data and uncertainty 
in the slope. 
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2. The origin of the intrinsic polarization 

Scattering and interstellar absorption cannot reproduce the observed wavelength dependence 
of the LP in Sco X-l, Cyg X-2 and GRO J1655-40, nor can they explain the short timescale vari- 
ability of the LP in Sco X-l. Polarized light is however expected from optically thin synchrotron 
emission from the jets in these systems, which has now been observed in a number of LMXBs in 
the NIR for both BHs [9, 22, 4, 35] and NSs [28, 33, 41]. The strength of LP is dependent on 
the ordering of the magnetic field structure near the emitting region, which for NIR optically thin 
synchrotron emission is close to the jet launch region at the base of the jet. Thermal emission from 
the disc dominates more at higher frequencies in the optical, whereas the jet is more dominant in 
the NIR, so the higher levels of LP in the NIR can be explained by a higher jet/disc flux ratio in this 
region of the spectrum. Moreover, the position angle (PA) of the polarization in Sco X-l (measured 
by us and by Shahbaz et al. 2008 [40]) is consistent with the magnetic field being aligned with the 
axis of the jet, which is known from radio observations [24]. We also considered LP caused by 
the scattering of the jets on the companion star surface, which can be ruled out because this would 
result in LP only at some phases of the binary orbit, which is not observed [32]. 

It is for these reasons that both us and Shahbaz et al. interpret this polarization as the first 
detections of the polarized inner regions of the jets. For a highly ordered magnetic field, the LP 
from optically thin synchrotron emission could be as high as ~ 70% [36, 2]. The levels detected 
(a few per cent at most) suggest a tangled magnetic field is present, unless the majority of the NIR 
emission is produced in the (unpolarized) disc or star. The disc and/or star do seem to contribute or 
even dominate the NIR for the three sources [29, 33], so /, the dimensionless parameter represent- 
ing the ordering of the magnetic field (LP = 0.7 x /), cannot be constrained for Sco X-l and Cyg 
X-2 because the jet/disc flux ratio is unknown. For GRO J 1655-40 we estimate the jet/disc ratio 
from broadband spectral energy distributions [29] and infer / = 0.41 ±0.19 [32]. To constrain the 
value of / further, we require data during times in which the jet is known to be dominating the 
emission. 

3. New NIR polarimetry of GX 339-4 

In 2008 September we observed the BH GX 339-4 on four dates (4 th , S th , II th and 25 th ) 
with ISAAC on the VLT (ESO Director's Discretionary Time programme 281.D-5029) in imaging 
polarimetry mode. At the time, the NIR and optical emission was highly variable (with high ampli- 
tude rms variability) and was correlated with the X-ray flux; the variability was found to originate 
in the jets [34, 6]. In Fig. 2 (left panel) we show the optical-NIR spectral energy distributions 
(SEDs) of GX 339-4 during this highly variable 'hard X-ray state' during which a steady compact 
jet is expected to exist (for state definitions related to jets see e.g. [16, 15, 1]). Data were acquired 
with ISAAC on the VLT, the Faulkes Telescope South (as part of a long-term monitoring campaign; 
Lewis et al., in preparation) and the South African Astronomical Observatory (SAAO) 0.75 metre 
telescope. The data were de-reddened assuming an extinction of Ay = 3.9 mag [23] and following 
the extinction law of [5]. 

Two components of emission are evident; one with a blue spectrum, a « +2.0 (where F v ©c 
v a ; possibly thermal emission from X-ray reprocessing on the accretion disc surface) and one 
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log (central frequency of filter; Hz) , Qg (central frequency of filter) 

Figure 2: Left panel: Optical-NIR spectral energy distributions of GX 339-4 during a low luminosity hard 
state in 2008 June - September. Data were acquired with ISAAC on the VLT, the Faulkes Telescope South 
and the SAAO 1.9m Telescope. The data were de-reddened assuming an extinction of Ay = 3.9 mag [23] 
and following the extinction law of [5]. Right panel: The rapidity of the variability from the Faulkes and 
VLT observations (the mean for each epoch in each filter), as a function of frequency of the filter (see text). 

with an approximately flat spectrum, a « 0.0 (from synchrotron emission from the jet; see also 
[9, 21, 6]). It is clear from these SEDs that the jet/disc flux ratio increases at lower frequencies. By 
extrapolating the average SED of the blue component into the NIR, we estimate that it contributes 
~ 24%, ~ 11% and ~ 6% of the average J, H and K s -band flux densities (however this is sensitive 
to the value of the extinction Ay, which is uncertain [3, 23]). By subtracting this contribution 
from the J, H and K s flux densities we estimate the jet to have mean flux densities of ~ 8.2 mJy 
in J, ~ 10.5 mJy in H and ~ 9.3 mJy in K s . The spectral index of the jet between J and H is 
a « —0.9, consistent with optically thin synchrotron, and between H and K s it is a « +0.4, an 
inverted spectrum inconsistent with optically thin synchrotron and more consistent with optically 
thick synchrotron (similar to one SED presented in [9]). 

In addition, the variability is strongest in the NIR J and H-bands (Fig. 2 right panel). The 
rapidity is defined as the mean of the differences between consecutive magnitudes, divided by the 
time resolution. The error bars represent the measured variability of a field star, which should be 
close to zero if noise is minimal. We do not include data in which the field star is measured to 
vary by > 1 milli-magnitudes per second. In the NIR J and H-bands, GX 339-4 varied by up to 7 
milli-magnitudes per second (up to 5 in K s ) whereas the maximum measured from the optical V 
and i' bands on six separate dates was < 3 milli-magnitudes per second (a similar trend was found 
by [17]). We expect the optically thin synchrotron emission to be highly variable and the optically 
thick emission less so, since the latter is a superposition of overlapping synchrotron components, 
which will dilute both the variability and polarization. Both the SED of the jet in the NIR and 
the wavelength-dependent variability are therefore consistent with the break (turnover) in the jet 
spectrum, between optically thick and thin emission, to reside around the H-band. This result is 
similar to other works of optical-NIR studies of GX 339-4 [9, 10]. 

In Fig. 3 we present the NIR polarimetric results. The Stokes parameters q and u 
(LP= \/q 2 + u 2 ) of GX 339-4 and a close field star are derived by correcting the measured values 
for instrumental polarization using the unpolarized standard star WD 1620-391. The LP is low for 
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-3 1 1 1 1 1 J ^ 1 1 1 1 1 1 1 

-3 -2 -10 12 3 -3 -2-10123 

Stokes q (%) Stokes q (%) 

Figure 3: Left panel: q and u Stokes parameters of GX 339-4, corrected for filter-dependent instrumental 
polarization (filled circles = K s -band, open diamonds = H, filled squares = J, open triangles = SZ). There 
is a correlation between q and u; the fit to the correlation (lcr errors) is shown ('NIR Pol. PA (fit)'), which 
corresponds to a PA of polarization of (17.0 ± 1.6)°. We also plot the correlation that corresponds to the 
measured PA of the resolved radio jet of (—64 ± 2)° [25] ('Radio jet PA). Right panel: q and u Stokes 
parameters of a close by field star, corrected for filter-dependent instrumental polarization. 



both GX 339-4 (LP < 2.8%) and the field star (LP < 1.3%). The values of q and u are clustered 
around zero approximately uniformly for the field star, which is expected for an ~ unpolarized 
source with some low level of noise. However for GX 339-4 there is a positive correlation between 
q and u. Using the Spearman's Rank method we find the correlation is significant at the 3.4<7 level 
and it also passes through (0,0); the fit and lcr error is u = (0.68 ±0.08) q +(0.002 ±0.001). This 
corresponds to a PA of polarization ( tan[2PA] = * ) of (17.0 ± 1.6)°. 

Since we know the jet dominates the NIR emission during these observations and the jet is 
variable on short timescales, we conclude that the polarization signature of the jet is detected, and 
changes rapidly. For optically thin synchrotron emission, the PA is parallel to the electric vector 
and perpendicular to the magnetic field vector. The measured polarization PA differs by (81 ± 4)° 
with the PA of the resolved radio jet of (—64 ± 2)° [25], so the magnetic field near the base of the 
jets in GX 339-4 is ~ parallel to the jet axis. 

On each date four measurements of LP were made (one in each of the four filters SZ, J, H 
and K s ); the time resolution is on the order of 3 - 4 minutes. If the emission is purely optically 
thin synchrotron, which it likely is in the H-J region of the SED (see above) then we estimate the 
ordering of the magnetic field to be / < 0.04 on these timescales. This is very low, and provides 
vital observational clues to the structure of the inner accretion (out)flows of BHs, which could be 
used to compare to those predicted by simulations and models of jet formation (these results will be 
published in Russell et al., in preparation). It is worth noting that it was predicted in 1982 [14] that 
the optical flares in GX 339-4 may be polarized; now almost 30 years later this has been confirmed. 
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4. New optical fast timing polarimetry of Sco X-l with simultaneous X-ray 

We observed the NS Sco X-l on two nights in February 2009 with the two channel High 
speed Photo-POlarimeter (HIPPO) [30, 31] mounted on the 1.9-m optical telescope of the South 
African Astronomical Observatory (SAAO). With HIPPO, Stokes parameters are recorded every 
0.1 sec and photometry every 1 millisecond, simultaneously in both channels. We used V and 
I-band filters for the two channels in order to achieve simultaneous flux and polarimetry at two 
wavelengths. 

X-ray observations of Sco X-l using the Proportional Counter Array (PCA) aboard the Rossi 
X-ray Timing Explorer (RXTE) satellite were acquired on seven consecutive dates, 19-25 Feb. 
On 19 th and 22 nd Feb, simultaneous optical-X-ray observations were arranged. This is the first 
time as far as we are aware, that optical polarimetry and X-ray observations of a LMXB have been 
performed simultaneously. Z-sources display a number of X-ray states which form a characteristic 
pattern in an X-ray colour-colour diagram or hardness-intensity diagram [19]. The states (the 
'horizontal branch', the 'normal branch' and the 'flaring branch') occupy different regions of these 
diagrams. It has been found that the behaviour of the radio jets correlates with the X-ray state. 
According to the unification of radio jets in NS Z-sources [20, 27], the jets are absent in the flaring 
branch, present in the horizontal branch and discrete ejections are launched in the normal branch. 

The X-ray and optical flux (and polarimetric) light curves are plotted in Fig. 4. The light 
curves from 19 th and 22 nd Feb are in the left and right panels, respectively. Comparing the two 
upper panels, the X-ray and optical fluxes were both brighter and more variable on 19 th compared 
to 22 nd . Integrating over 25 minutes, the LP is measured to be < 1 % (also shown in the upper 
panels) but significant, and consistent with an interstellar origin (see Section 2 and Fig. 1). In the 
centre two panels of Fig. 4 we plot measurements of polarization but binned over only 30 seconds. 
With smaller bins the error bars are larger, but we plot just the 3<7 detections of polarization in these 
centre panels. We expect some false 3a detections; below we discuss whether (some of) these 
detections are real or are due to noise. LP is again plotted in the lower panels, this time from just 
10-second bins and we only plot the data that are 3d detections in both channels simultaneously. 
In the lower right panel (22 nd Feb), there appears to be a polarization flare lasting 2 bins (~ 20 sec). 

In Fig. 5 the X-ray hardness-intensity diagrams are plotted. We compared the data on Feb 
19 th (black data in left panel) and Feb 22 nd (black data in right panel) with those of the period 
Feb 19-25 r/z , and compared the X-ray power spectra with those in the literature [42]. We find the 
source was in the 'flaring branch' on 19 th and the 'normal branch' on Feb 22 nd of the NS Z-source 
classifications. 

To test whether these optical polarization flares could be real, we plot the PA of the significant 
(at the 3<7 level) detections of LP (in red) and the PA of the non-detections (black data) in Fig. 6 
(using 10-sec binned data). All non-detections are ~ evenly distributed in PA. On Feb 19 th (left 
panels), almost all 3<7 detections of LP appear to have random PAs (again approximately evenly 
distributed) but on Feb 22 nd (right panels) most of the 3 a detections are clustered around PA ~ 
100°. When we detected polarized NIR flares from the jet in 2006 (section 2), their PA was also 
~ 100° (shown as the grey regions in Fig. 6), which suggests the flares of LP could be real, and 
originate in the jets. 

We see the polarized flares during the normal branch and not (significantly) during the flaring 
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Figure 4: Simultaneous light curves: optical flux and polarization in two channels (V-band in blue and 
I-band in red; from HIPPO) and X-ray 3-60 keV flux per PCU (in black; RXTE PCA). Left panels: Data 
taken on 2009-02-19. The time resolution of the optical flux and X-ray flux is 5 seconds and 1 second, 
respectively. The polarization measurements are made taking data over timescales of 25 minutes ('25 min 
bins'; top panel), 30 seconds ('30s bins'; only 3d detections of polarization in one of the two channels 
are shown; centre panel) and 10 seconds ('10s bins'; only 3<7 detections of polarization in both channels 
simultaneously are shown; lower panel). The X-ray light curve and '25 min bins' optical polarization data 
are also shown in the centre and lower panels for comparison. Right panels: Same as the left panels but for 
data taken on 2009-02-22. 
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Figure 5: X-ray hardness-intensity diagram for Sco X-l. The count rate is in the 3-60 keV range and the 
soft colour is defined by the 4.4-6.2 keV flux divided by the 3.0-4.4 keV flux. All data from the seven 
dates are shown in grey, and the data on 19 th and 22 nd Feb are shown in black in the left and right panels, 
respectively. 
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branch; this is consistent with the global picture of jet-disc coupling in NS Z-sources [27]. The 
implication is that the optical emission from Sco X-l is sometimes polarized by several per cent, 
but on timescales of 10 seconds or less. The jets appear to cause this transient polarization, and 
the flares could arise due to brief ordering of the predominantly tangled magnetic field near the jet 
base, or intermittency of the jets themselves. 

5. Conclusions 

The polarimetric signature of synchrotron-emitting jets can be detected in the infrared and op- 
tical in LMXBs. The results so far (from at least one BH source, GX 339-4 and one NS source, Sco 
X-l) suggest that a generally turbulent (only partially ordered), rapidly changing magnetic field is 
present near the base of the jets close to where they are launched. When significant polarization 
is measured, the position angle is consistent with a magnetic field aligned with (parallel to) the jet 
axis. The polarimetric variability timescales are at least as short as 10 seconds in Sco X-l, and the 
net polarization signal is diluted in exposures longer than this timescale. The short timescale of 
the flares of polarization places observational constraints on models for the jet launching process 
and the local magnetic field configuration. The magnetic field configuration, and hence polariza- 
tion, cannot vary on timescales shorter than the light travel time across the emitting region. It is 
therefore very likely that a HTRA instrument with polarimetric capabilities on an extremely large 
telescope such as the planned European Extremely Large Telescope (E-ELT) will uncover the flux 
and polarization properties of these objects over all plausible variability timescales. 
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